How many electrons are needed to flip a local spin? 
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Considering the spin of a local magnetic atom as a quantum mechanical operator, we illustrate 
the dynamics of a local spin interacting with a ballistic electron represented by a wave packet. This 
approach improves the semi-classical approximation and provides a complete quantum mechanical 
understanding for spin transfer phenomena. Sending spin-polarized electrons towards a local mag- 
netic atom one after another, we estimate the minimum number of electrons needed to flip a local 
spin. 
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Electrons generally interact with magnetic atoms 
through a spin-flip interaction. For example, this is a 
fundamental mechanism for spin transfer 0,13] from spin- 
polarized electrons to a magnetic moment in a ferromag- 
net Spin transfer results in a classical torque 

exerted on the magnetic moment and, thus, enables us 
to control a local structure in a ferromagnetic film. This 
picture is semi-classical and seems to work well n a practi- 
cal sense 0,0, However, this view also causes some 
conceptual difficulties, and cannot answer fundamental 
questions associated with spin transfer. For example, it 
seems intuitive that the stronger the spin-flip interaction 
is, the easier it is to flip a local spin through the inter- 
action. However, a more rigorous quantum mechanical 
treatment of the problem will illustrate that this intu- 
ition is incorrect. 

The goals in this paper are following: (1) we scrutinize 
the physics of spin transfer between an incoming elec- 
tron and a localized magnetic atom with spin S, i.e. we 
calculate the expectation value of the local spin opera- 
tor as a function of time (t) and illustrate the physics in 
detail. (2) We provide an estimation of how many elec- 
trons are needed to flip the local spin by sending spin- 
polarized electrons one after another with an interval suf- 
ficiently long that no interference between two electrons 
occurs. For illustrative purposes, we consider spin values 
of S = 1/2, 1, and 3/2 to represent the local spin. 

Imagine a magnetic atom residing at the origin of the 
X axis with spin pointing to the negative Z axis. Its spin 
state can be represented by | — S) at t — 0. We consider a 
normalized wave packet <j> of an electron with spin along 
the positive Z direction away from the origin. Then, the 
total wave function is initially 'J? = <f>\+)\ — S), where 
|+) is the spin state of an incoming electron. Since the 
local atom is neutral, the Coulomb interaction will not be 
included. This problem can be set up in one dimension. 



The Hamiltonian we consider is 



H = — 2 Jq<t ■ SS(x) 

2m 



(1) 



where p 2 /2m is the kinetic energy of the incoming elec- 
tron, a is the electron spin operator, and Jo (> 0) is the 
coupling of the spin-flip interaction. We use units such 
that h = c = 1. Let us introduce A/a = 2m Jo, where a 
is a typical length scale in the problem and will be set to 
unity. The initial wave packet is normalized and given 
bycf>(x) = (2 7ra )- 1 /4 e4 fco(^ ) e -(x+ a; o) 2 /4a, Such a wave 

packet describes an "electron" with mean position — xq 
and mean momentum kg. The uncertainties associated 
with the packet are Aa:o = %/2a and Afco = 1/V2a. In 
order to calculate the expectation value of the Z com- 
ponent of the local spin (S z (t)) = ($>(x, t)\S z \&(x, t)), 
where (• • •) means an integration over x, we need to know 
the time evolution of the total wave function ^(x,t). 
Since the initial state ^(x, 0) is not an eigenstate of the 
Hamiltonian because of the spin-flip coupling, ^(x,0) 
should be decomposed into two channels depending on 
the total spin j± = S ± 1/2. 

The Schrodinger equation to solve is 



c^{x) + 2mEV(x) = 



-SX5(x)V(x) for j + 
(S + l)X5(x)^(x) for j_ 

(2) 

A wave function with j + sees a potential well while a 
wave function with j_ feels a potential barrier. The time 
evolution of each channel is different because the eigen- 
states of each equation are different, and the overlap of 
the two wave functions determines the dynamics of the 
local spin as we will soon show. The time evolution of 
the wave function is a combination of the time evolutions 
of the two channels of j± : 



9(x,t) = <t>+(x,t) 
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\j+,-S + 1/2) + 4>-(x,t) 



2S 
2S+1 



\j-,-S +1/2) 



(3) 
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where 4>±{x,t) is the time evolution of the wave packet 
in the presence of the well/barrier, and </>±(a;,0) = </>(x). 
The Z component of the total spin is conserved to be 
— 5 + 1/2. The expectation value of the Z component of 
the local spin as a function of time is calculated to be 



(SM) = 



s 



(2s + iy 



{45 2 + 45 + 4Re[(0+|0_)] -3} 



(4) 

Using the Schwartz inequality: |(0 + |0_)| 2 < 1, we find 
a constraint of (S z (t)) as follows: —5 < (S z (t)) < 
-S [45 2 + 45 - 7] /(25 + l) 2 . This constraint is more 
restrictive than the obvious one that requires a max- 
imum change in the spin to be unity since the spin 
of the incoming electron is 1/2. Note that there al- 
ways exists a probability that the state of the local 
spin remains unchanged even after the interaction; there- 
fore, actual values of (S z (t)) should be evaluated quan- 
tum mechanically. We can also determine the local 
spin state from the total wave function ^f(t). The 
spin state is written as oi| — 5) + b\\ — S + 1), where 
|ai| 2 = {l + 45 2 + 45Re[(0 + |0_)]}/(25 + l) 2 and 
|6i| 2 = 45 {1 - Re /(25 + l) 2 . It is straight- 

forward to show that (S z (t)) = — 5 + |&i| 2 . Without loss 
of generality we can assume that a\ and b\ are real. 

As mentioned for (S z (t)), we need to solve Eq. J5J) be- 
cause the time evolution of the wave packet <fi± will be 
governed by the solutions of this equation. Let us con- 
sider a Hamiltonian H± = p 2 /2m + (X±/2m)5(x), where 
A+ = — 5A for the j + channel while A_ = (5 + 1)A for 
the other channel. This Hamiltonian could be treated 
as a scattering problem for E > to obtain an asymp- 
totic solution which is plane wave-like. Note that the 
asymptotic solution does not describe the detailed dy- 
namics of the local spin. However, it implies important 
physics associated with the problem. The asymptotic so- 
lution is r) k (x) = [e lkx + R±e- tkx ] G(-x) +T±e ikx Q(x), 
where 0(x) is a step function. The momentum is well- 
defined while the position cannot be; namely, Afc = 
but Ax = oo. The transmittance and the reflectance 
are determined by the boundary condition at x = 0, 
and they are |T±| 2 = 4/ [4 + (A±/fc) 2 ] and |i?±| 2 = 
(A±/fc) 2 / [4 + (A±/fc) 2 ] . It is worth mentioning that 
|T±| 2 and |i?±| 2 depend only on (A±/fc) 2 ; furthermore, 
for large A±/fc both channels give unit reflectance and 
almost zero transmission. As we showed in Eq. (0J, 
the dynamics of the local spin depends on the overlap 
of the two wave functions 4>+{x,t) and </>_(x, i). The 
above exercise indicates that for a large A, <j>± would 
not differ significantly from each other. In this instance, 
Re [(</>+ 1 (/>_)] « 1; consequently, (5 2 ) w —5. In other 
words, if the spin-flip coupling is very large, it becomes 
more difficult to flip a local spin. This seems to oppose 
the semi-classical understanding but we will show that 
this is the case. We should mention that the X±/k scal- 
ing in |T±| 2 and |i?±| 2 takes place because of zero un- 



certainty in the momentum (Afc = 0) for a plane wave. 
Rigorously speaking, however, one cannot expect such a 
perfect scaling when a wave packet with Afc ^ is used 
instead of a plane wave. 

We obtain the eigenstates and the corresponding eigen- 
values of H± introducing a box of the length 2L (-L < 
x < L) with a periodic boundary condition: tp(—L) = 
tp(L) and d x ip(—L) = d x tp(L), where tp(x) is an eigen- 
state of H± with an eigenvalue E±. Since the potential 
is symmetric about x = 0, the eigenstates are either even 
or odd; ip e ^ ± (x) or ip , p± (x). We found 0] for the even 
solution 



i>e,k ± (x) 



1 



cos(fc±x) + — -^sin(fc±|x| 
2fc± 



(5) 



where N k± = L [l + (A±/2fc±) 2 ] + A±/2fc|, and for the 
odd solutions ip 0tP± (x) = sin(p±x)/yL. The correspond- 
ing eigenvalues are Ek ± — k 2 _/2m, where k± is a solu- 
tion of tan(fc±L) = A±/2fc±, while E p± = p±/2m with 
p± = n±ix J L (n± is an integer). Since the potentials 
do not affect the odd solutions, the dynamics of the lo- 
cal spin is determined only by the even solutions. For 
A + (< 0), the number of the even states decreases by 
one and a single bound state occurs for E < to keep 
the total number of eigenstates unchanged. As long as 
we initially put a wave packet far away from x = 0, the 
bound state does not participate in the time evolution of 
the wave packet ^l|- Now the time evolution of <f>±(x, t) 
for both channels can be written as 

4>±(x,t) = ^e-^± t C fe± ^ e , fc± (x)+^e-^± t C p± V ,p ± 

k± p± 

(6) 

where C k± = (ip e ,k±\4>) an d C p± = (ip 0!P± \(j)). Using 
these expressions we obtain the overlap between <p + and 
6- as follows: 



.|*_) = (25 + l)A]T 



e i(E h+ -E k _)t ^ ) { ^ 



k 2 



l r 

2 



l 



-(fco/Afco) 2 -(*o/Azo) S 



(7) 



where Afco and Axo are the uncertainties associated with 
the initial wave packet. Note that the time dependence 
of (0_|_|^>_) is determined only by the even solutions while 
the odd solutions contribute to a constant, which would 
be close to 1/2 if the momentum or the position is well de- 
fined and finite initially. Since we consider a wave packet 
far away from the origin with a finite mean position at 
t = 0, the second term in Eq. (JJJ is approximately 1/2 
to high precision. 

Fig. 1 illustrates the time evolution of the wave packet 
for the j_ channel away from a local spin of 5 = 1/2. 
We introduce a dimensionless time r = t/2ma 2 . Initially 
the wave packet is located at x = — xo (xo = 100) with 
a mean momentum fco = 1. We choose a to be 10 and L 



3 



to be 10 3 ~ 10 4 . As long as L ^> xq, L is not important. 
The wave packet does not interact with the potential with 
A_ = (5 + 1)A = 3.45 until r ~ 30. At r = 50, the po- 
tential strongly scatters the wave packet and a part of 
the packet is transmitted to the other side of the poten- 
tial. After t = 150, the reflected and transmitted parts 
of the packets freely move away from the potential in the 
opposite direction as indicated by the arrows. The time 
evolution for j + is similar and not shown in this figure. 
In the insert, we show (S z (t)} as a function of r. It in- 
creases from (5 2 (0)) = -0.5 to (5 Z (150)) ~ -0.01, and 
then becomes saturated. The saturation is expected be- 
cause after r = 150 the wave packet does not interact 
with the potential any more. In fact, the mean momen- 
tum determines when (5 2 ) becomes saturated. We define 
the saturation value of (5 Z ) as Si, which depends on kg 
and A. It is this value that we use to estimate how many 
spin-polarized electrons are needed to flip a local spin as 
we will show later. 

It is possible to control the mean momentum reason- 
ably well while the spin-flip coupling is uncontrollable 
and even not well known experimentally. In this sense, it 
does not seem possible to pinpoint how many electrons 
are needed to flip a specific local spin. Nevertheless, it is 
theoretically feasible to estimate the minimum number of 
electrons to flip a local spin by calculating the maximum 
value of Si denoted as S™ ax . To this end, we evaluate Si 
as a function of kg and A because a scaling relation is not 
expected to hold for the wave packet with the momentum 
uncertainty (l/y/2a). We choose the same value of a as 
before. Fig. 2 is a contour plot of Si in the (A, ho) plane. 
Interestingly, we found a X/ko scaling holds fairly well in 
this case as well. Moreover, 5™ a;r of a local spin S = 1/2 
is approximately zero for A ~ 2.3/co; for example, A = 6.9 
with feo = 3. We obtained the same scaling behavior 
for other values of a. Note that this value of S™ ax is 
smaller than the mathematical upper bound based on the 
Schwartz inequality, which is a half for S = 1/2. Looking 
in particular at the lower part of Fig. 2, we know that for 
a given kg, Si decreases with increasing A. This clearly 
indicates that if the coupling is too large, it becomes more 
and more difficult to flip the local spin, which is consis- 
tent with our early analysis. We plot Si for the local 
spin 5=1 and S = 3/2 in Fig. 3 as a function of X/kg 
with a — 10. There is also a similar scaling behavior in 
Si while different values of X/kg give S™ ax . For 5=1, 
the maximum value is S™ ax ~ —0.556 (< —1/9) along 
A ~ 1.4fc , and 5]™" - -1.126 (< -3/4) along A ~ k 
for S = 3/2, where —1/9 and —3/4 are the mathematical 
upper bounds for 5=1 and 5 = 3/2, respectively. 

Now consider the spin-polarized electrons sent towards 
the local spin one after another over an interval To which 
is sufficiently long so as to prevent any interference be- 
tween two electrons. That is, we wait long enough before 
sending the second electron so that the first has cleared 
away from the region of interest, and left the local spin in 



astateai|— 5)+6i|— 5+1), where a i and&i have acquired 
saturated values after a time r « 200 (see Fig. 1). Note 
that this time can be shorter if we increase the mean mo- 
mentum, for example. At t = t , the total wave function 

is ct>\+)[ai(T )\-S) + bi(r )\-S + l)]. Since |+)|-5) 

and |+) | — 5 + 1) are not eigenstates of the Hamiltonian 
in general, we need to decompose each state into two 
channels for j± with appropriate Clebsch-Gordan coeffi- 
cients ^3 • Later the spin state of the local spin will be 
a 2 \- 5) + b 2 \- 5 +1) + c 2 |- 5 + 2) . The coefficients a 2 , b 2 , 
and c 2 can be determined from the total wave function, 
and they turn out to be functions of ai and b%. Repeat- 
ing this procedure we can calculate the expectation value 
of the local spin after sending the n-th electron. 

Using 5 = 1/2, 1 and 3/2 we illustrate the procedure 
and estimate the minimum number of electrons to flip a 
local spin. The expansion for an arbitrary spin can be 
done systematically. For a local spin 5 = 1/2, j + = 1 and 
j_ = 0. After the first electron interacts with the local 
spin, the total wave function is <3> = -y=<^i|10) + ^<^o|00). 
The spin state of the atom is <ii| — 1/2) + &i|l/2), where 
a\ = l/2-5i and b\ = l/2+5 x , and (S z ) = Si. When we 



send the second electron, ^ 
Mi|H) 

where a| = a\ and b\ = (a\ 



.ai0i|lO) + 4-ai<fo|OO)- 



The local state becomes a 2 \ — 1/2) + 62 1 1/2), 
= a\ and b\ = (a\ + l) b\, and 5 2 = 1/4 + Si - 
Sf, where S n is (S z ) after the n-th electron. The third 
electron gives 5 3 = 1/4 + (Si + 5 2 ) /2 - 5i5 2 , and the 
n-th electron leaves S n = 1/4+ (Si + 5„_i) /2 — SiS n -i- 
We can express S„ in terms of Si as follows: 



Sn -2~\2~ Sl 



(8) 



Using a similar procedure we obtain 5„ for 5=1 and 
5 = 3/2 

S n = l-?i(-5i)" _1 + (?i-2)(-5i)" for 5 = 1 



5„ — - - (-Si 



1/2)™ (3 + 4n)5i +3(1 + An)/2 



6(-45i/3-lf for 5 = 3/2 



(9) 



Note that 5o = —5, which means that initially the lo- 
cal spin state is | — 5) while S„ — > 5 as n increases 
for a given Si; in other words, the local spin becomes 
flipped. Mathematically speaking, 5 n = 5 only when 
n = 00. This is because there always exists a quan- 
tum mechanical possibility that the spin state remains 
unchanged even after the spin-flip interaction. Never- 
theless, when S n becomes sufficiently close to 5, we can 
claim that the local spin has flipped. When Si = S™ ax , 
we can approximate S™ ax ~ 5 — 25 e~° sn , where 
Pi 12 — 0.7, Pi ~ 0.35, and /3 3 / 2 ~ 0.23. Let us define 
the minimum number N$ to be the number which satis- 
fies, for example, 25 e~@ sNs = 10~ 5 . Then we evaluate 
N s = [51n(10) + ln(25)] //3 S - Our estimation shows that 
the minimum number of the spin-polarized electrons to 
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flip a local spin of 1/2 is about 16; namely, JV1/2 — 16. 
For S = 1, Ni ~ 35 while -/V 3 / 2 ~ 55 for the local spin of 
3/2. Therefore, less electrons are needed to flip a smaller 
spin as one might expect. 

In summary, using straightforward quantum mechan- 
ics, we have studied the time evolution of the spin of a 
local magnetic atom under a spin-flip interaction with 
an incoming electron. This treatment goes beyond the 
semi-classical approximation, which considers the local 
moment as a classical vector. The expectation value of 
the spin operator has been evaluated using the wave func- 
tion of the electron, which is the solution of the time 
dependent Schrodinger equation. Sending spin-polarized 
electrons towards a local magnetic atom one after an- 
other, we also provide an estimate of how many elec- 
trons are needed to flip a local spin. For an experimen- 
tal realization of our estimate, we suggest a setup where 
a magnetic atom is fixed at the hub of a wheel, while 
spin-polarized electrons are sent towards the atom along 
orthogonal spokes in the wheel. 
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Fig. 1 ( Kim et al ) 
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FIG. 1: (Color online) The time evolution of the wave packet 
for the j- channel. The other channel shows similar behavior 
and is not plotted. The local spin is 1/2. The inset describes 
dynamics of the local spin. After r = 150, (S z ) becomes 
saturated. 
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FIG. 2: (Color online) Contour plot of Si for S = 1/2 in 
the (A, ko) plane. Note a (A/feo) scaling behavior with a 
maximum value of Si along A = 2.3fco. 



Fig. 3 ( Kim et al ) 
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FIG. 3: (Color online) Si as a function of X/ko for 5 = 1 
and S = 3/2, for various values of fco- The maximum value 
of Si is about -0.556 and -1.126 for S = 1 and S = 3/2, 
respectively. 



